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Abstract 

Observation indicates that the expansion of the Universe is accelerating and 
favours a dynamical cosmological constant, A(t). We consider the possibility 
that this is due to a scalar field which has undergone a very recent phase transi- 
tion. We study a simple class of model, corresponding to a cp^ potential with a 
time-dependent mass squared term. For the models considered the phase transi- 
tion occurs at a red shift z < 1.2. The evolution of the equation of state uj^ and 
energy density with time is distinct from existing dynamical A models based 
on slowly rolling fields, with uj^ and p,^ rapidly changing in a characteristic way 
following the transition. The (p energy density is composed of a time-dependent 
vacuum energy and coherently oscillating condensate component with a nega- 
tive pressure. The condensate component will typically collapse to form non- 
topological soliton lumps, '(/)-axitons', which smoothly populate the Universe. 
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Perhaps the most remarkable cosmological observation of recent times is the accel- 
erating expansion of the Universe . This requires the existence of an energy density 
which has a negative pressure and which is smooth on the lOMpc scales relevant to 
dynamical estimates of the density of conventional dark matter. Other evidence for 
this comes from estimates of the age of the Universe, the Hubble constant, the baryon 
fraction in clusters, the galactic power spectrum and CMB measurements 0, ^. The 
simplest explaination for a smooth energy density with a negative pressure is a time- 
independent cosmological constant, A. However, there are a number of problems with 
a fixed cosmological constant. The first is that it is difficult to understand why the 
cosmological constant should just be dominating the energy density at the present 
epoch [^. This has led to two different views. One is that the cosmological constant is 
due to the evolution of a scalar field ('tracking solution' see also P) whose energy 
density becomes significant recently due to dynamical effects. In general this approach 
has difficulties with nucleosynthesis and the present equation of state (A particu- 
larly promising version which may overcome these problems is given in [||.) The other 
view is that the cosmological constant has become dominant in the present epoch due 
to anthropic selection (AS) 0]. However, even if we assume that AS is responsible for 
the dominance of the cosmological constant, there are still problems. It is difficult to 
understand why a fixed A is so small compared to the mass scales of particle physics 
(A ~ 10~^^°Mpi). This has led to the suggestion that the absolute minimum of the 
vacuum energy should be exactly zero and that the smooth energy density is due to 
the evolution of a scalar field towards the minimum, resulting in a dynamical cosmo- 
logical constant, A(t) |10|, |T3|, 0. (An interesting alternatve based on quantum 
spinodial fiuctuations has been suggested in |jl5|.) There is also observational evidence 
in support of a specifically dynamical cosmological constant [|1^. The amplitude of 
the COBE-normalized galaxy clustering power spectrum is too large in the case of a 
fixed A. However, if the effective cosmological constant is decreasing with time, then 
the amplitude of the power spectrum is reduced, alleviating the problem ||12|. 

There have been a number of suggestions regarding the nature of the dynamical 
cosmological constant. Most popular are models based on pseudo-Nambu-Goldstone 



1 



bosons (PNGBs) [jlT], |T^, exponential potentials ||, |T3[ and inverse-power 
law potentials |]14|, |10]. All of these produce characteristic time-dependent energy 
densities and pressures (= oj^p^) which come to dominate at recent red-shifts. 
For the PNGB models, the pressure slowly tends from negative values towards zero. 
For tracking models based on exponential and inverse-power law potentials, on the 
other hand, the pressure is evolving from a value typically positive and close to zero 
(corresponding to matter tracking) towards a negative value. This should allow the 
models to be distinguished by precision CMB angular spectra measurements |13|. All 
of these models are based on very light scalar fields [m^ ^ iJo, where Ho is the present 
value of the expansion rate) which are slowly rolling at the present time, |0/0| ^ Ho- 

In the present paper we wish to introduce an alternative model for a dynamical 
cosmological constant. This is based on the idea that the negative pressure energy 
density is associated with a conventional metastable false vacuum. This is perhaps the 
simplest form of negative pressure energy density in particle physics models. However, 
in order to have a dynamical cosmological constant we require that a phase transition 
from the metastable phase has occured recently. To achieve this we will consider the 
scalar field (f) to have a time dependent mass squared term which has recently become 
negative. We will refer to this scheme as the 'very recent phase transition' (VRPT) 
scenario for a dynamical cosmological constant. The VRPT scenario is quite distinct 
from previous models based on slowly rolling scalar fields. In particular, there will be 
no need for the mass scale of the scalar to be extremely small, and the pressure and 
energy density will be rapidly evolving in a characteristic way at recent times following 
the phase transition. This should allow the VRPT scenario to be distinguished from 
the others by precision CMB measurements. 

We will consider the usual spontaneous symmetry breaking potential for a real 
scalar field with — symmetry^, 

n0) = -^0^ + ^0^ + A ; A = ^, (1) 



^This can be generalized to a complex scalar. For a real scalar field domain walls may form during 
the VRPT, whilst for a complex field global strings may form. 
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where 

and Qc is the scale factor at the time of the transition. In general the VRPT scenario 
requires an additional time-dependent cj? term, which we will refer to as a 'stabilizing 
interaction'. We will discuss some possible sources for the stabilizing interaction later, 
but for now we simply model it phenomenologically. We will see that such time- 
dependent mass squared terms with integer n can arise naturally in plausible models. 

Typically the mass scale of the potential will be very large compared with Ho- 
Thus following the phase transition (f) will be coherently oscillating about the relatively 
slowly evolving time- dependent minimum of its potential. In order to discuss the time 
evolution of the energy density and equation of state, we follow the discussion of 



1^, based on averaging over the rapid oscillations of the scalar about the minimum 
of its potential. 

The (p equation of motion is given by 

dV 

<p + ?>H<j) = -— . (3) 



This may be rewritten as 

where the partial derivative is with respect to constant yU^(t).0 Taking the time average 
over an oscillation cycle, we obtain 

at 

where 



-3H-fp^ , (5) 



^ Ul - V/V^a.)-^'^d<P • ^ > 

denotes integration over one oscillation cycle and Vmax = P4> is the maximum 
energy density during an oscillation cycle. 7 corresponds to the time average of 0^ 



^In practice we evolve the energy density and pressure of 4> by incrementing the energy density 
with held constant and then incrementing and calculating the pressure. Therefore we do 
not use the total derivative with respect to t given in |16| for an explicitly time-dependent potential. 
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over an oscillation cycle. The equation of state of the energy density is then given 

by 

uj^=p^/p^ = -i -I . (7) 



We have calculated the time evolution of the energy density and equation of state for 
different values of n as a function of red-shift, z. The parameters of the potential are 
chosen such that we have at present ~ 0.7 and Vim ~ 0.3 (where Vim is the density 
of conventional clustered dark matter). 

The numerical solutions for and as a function of z depend on pi and A 
only through the ratio p'^/X'-, we have calculated the evolution for the case A = 1. 
Because of this scaling property, for a given value of today, uj^ o, the evolution of 
the (f) equation of state and energy density is completely fixed by the value of n in the 
stabilizing interaction. 

In Figure 1 we show the evolution of the equation of state uj^ as a function of 
red-shift for n = 2, 3, 4. In Table 1 we give the parameters of the VRPT {oj^ o, Po (for 
A = 1) and Zc, the red-shift at which the transition occurs) for the cases n = 2 and 
n = 4. For n > 2 and cu,/, o ^ —0.6 we find that the phase transition occurs at z^ < 1.2. 
(We consider uj^o ^ —0.6, in keeping with observational limits for the case of a fixed 
LOfj, Since in our case u^p is decreasing with z, this limit should be conservative.) 
From Table 1 we see that the mass of the scalar is ^ 10^^ eV for A 1. In Figure 2 
we show the evolution of the energy density together with the matter energy density 
flm for the cases n = 2 and n = 4, where we have normalized the energy density by 
taking the ratio to the present critical density pc- Both cu,^ and p,^ rapidly change at 
recent red-shifts, the more so for larger values of n. 

In addition, we have considered the effect of the VRPT on the age of the Universe, 
given by 

tu = \H~'fu ; fu = ll^ + ^m{t)V . (8) 

The age of globular clusters requires that fu = 1.5 ±0.3 [@, 0. For a fixed cosmological 
constant and = 0.7, fu = 1.45. For the VRPT fu is generally smaller, but not 
significantly so. The largest deviation in the examples considered corresponds to n = 2 
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Figure 1: Evolution of the equation of state for different n with uj^ = —0.6 and 
uj^ = —0.8 at present. 

and uj(i,o = —0.6, for which fu = 1.38. 



Table 1. VRPT parameters. 



n 


^(j, o 






2 


-0.6 
-0.8 


3.45 X 10"^ eV 
3.15 X 10-3 eV 


1.20 
0.47 


3 


-0.6 
-0.8 


3.27 X IQ-'^ eV 
3.13 X 10-3 eV 


0.57 
0.28 


4 


-0.6 
-0.8 


3.21 X 10^3 gv 
3.10 X 10^3 gv 


0.38 
0.20 



After the phase transition has occuredQ, the energy density will be composed of a 
time dependent vacuum energy p^jac and a (p condensate component posc corresponding 
to coherent oscillations about the time dependent minimum of the potential. In Figure 

^We are assuming that (j> is out of thermal equilibrium, so that the transition is purely dynamical 
in nature. 
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Figure 2: Evolution of the energy density Q^(t) = p^{t)/pc for ri = 2 and n = 4. 



3 we show the equation of state and energy density of the (p oscillations as a function 
of z. We see that the equation of state rapidly tends towards the value Uosc = corre- 
sponding to effectively 0^ oscillations about the minimum. Nevertheless, the equation 
of state and so the pressure is significantly negative throughout. In Figure 4 we show 
the equation of state and energy density associated with the time dependent vacuum 
energy, p^ac = ^(0mm(^))- We see that the equation of state does not significantly 
deviate from the value expected for a constant vacuum energy density, ua = —1. 

The negative pressure of the oscillations implies that the condensate is unstable 
with respect to spatial perturbations of |T^. Such perturbations may be expected 
to exist, coming, for example, from thermal or inflationary quantum fluctuations. In 
the linear approximation the spatial perturbations evolve as 



(9) 



(where cu^ < is assumed). Therefore, 

50k = exp(|u;0|^/^|k|t)(50ko 



(10) 
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Figure 3: Evolution of the equation of state and energy density of the coherent 
oscillations. 
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Figure 4: Evolution of the equation of state and energy density of the time dependent 



vacuum energy. 
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This is true so long as the wavenumber of the perturbation satisfies |kp ^ 41^7^17710 [p!9[] , 
otherwise the positive pressure associated with gradient energy in the perturbation will 
overcome the negative pressure responsible for the growth of the perturbation. (For 
this reason, negative pressure effects play no role in the dynamics of A(t) models based 
on slowly rolling fields, since their small mass, m<^ < i^o, implies that perturbations 
on sub-horizon scales cannot grow.) The perturbations in the condensate will grow 
exponentially until they become non-linear. The condensate will then fragment into 
non-topological soliton lumps which we will refer to as '0-axitons' (following the ex- 
istence of similar objects in axion cosmology |^ and Affleck-Dine baryogenesis ||21|). 
The radius of the 0-axitons will be determined by the first perturbation mode to go 
non-linear, «i (Icj^l^/^m^)"^ (assuming the exponential factor to be dominant in de- 
termining non-linearity). This will occur in a time 5t ^ (|u;(^|m0)~^log((/)/50k) ^ H^^. 

Thus shortly after the VRPT the Universe will typically be filled with 0-axitons. 
So long as -C lOMpc, which is true if ^ 10~^^ GeV, the 0-axiton density will 
initially act as a smooth component of pressureless matter as far as determinations of 
the dark matter density are concerned. 

We need to check that subsequent infall into galactic halos does not result in the 
0-axiton density clustering on lOMpc scales and so no longer being effectively smooth. 
We do this via a simple Newtonian argument. We consider the mean distance between 
galaxies to be Rgai ~ lOMpc, with 0-axitons being smoothly distributed initially. The 
time scale for the 0-axitons to fall a distance ~ Rgai due to the attraction of a galaxy 
of mass Mgai is then 

^3/2 

tinfall ~ , • fill 



Using Mgai ~ IO^^Mq for the mean galactic mass, this gives ti„/a// ~ 1-6 x lO^^yr, 
which is much longer than the age of the Universe, tu ^ 10^'^yr. Thus infall will not 
significantly alter the smoothness of the 0-axiton distribution. (The same argument 
holds for the (p particles in a homogeneous condensate.) Although we have estimated 
this at the present time, it holds for earlier times also, since tu (oc H'^) and tinjaii (oc 

-,3/2- 



^lai) both proportional to c?^'^ . 
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Figure 5: Evolution of the total equation of state and energy density for the case of 
pressureless 0-axitons. 

Therefore in the presence of spatial perturbations, the energy being fed into the 
coherent oscillations will be converted into a smooth density of pressureless 0-axitons. 
This will slightly alter the evolution of the total equation of state and energy density 
from the case where there are coherent oscillations with negative pressure. In Figure 
5 we show the effect of replacing the condensate with pressureless 0-axiton matter. 
The effect is a small alteration of the total equation of state and energy density as 
a function of z. 

After the phase transition, the energy density consists of a time-dependent vac- 
uum energy 



4A 



(12) 



with cuvac = — 1 to a very good approximation, and an energy density either in the form 
of a coherently oscillating field with a negative pressure or, more likely, in the form 
of pressureless 0-axitons. However, the total uj^ is insensitive to the pressure in the 
condensate component, as seen from Figure 5. This means that a;^ > — 1 following the 
VRPT is simply due to the dilution of cu^ by the low pressure component of smooth 
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(j) condensate or 0-axiton matter, 

^ oscPosc ~l~ ^vacPvac Pvac o\ 

^4> = 7 ~ 7 • (13) 

Pose ~r Pvac Pose ~r Pvac 

So far we have introduced a simple time-dependent mass squared term ('stabilizing 
interaction') in order to trigger the phase transition and so produce a dynamical cos- 
mological constant. We now consider some ways in which this term could be generated 
in particle physics models. One possibility is that the field might couple to a light 
field which has a thermal distribution, such as the photons or neutrinos. For example, 
one could consider a coupling of the form -^Cke, where Cke = —\F^''F^y, ipd^'-j^il). 
On taking the average over thermal fiuctuations of the fields in Cke we obtain an ef- 
fective cjP' term ~ -^^0^. With T oc CL ^ this results in a stabilizing interaction with 
n = 4. Alternatively we might consider a very light additional scalar field x with a 
thermal distribution (and a small enough energy density so as to avoid nucleosynthesis 
constraints), coupling to via a term Oii averaging over thermal fiuctuations 

of the X field, this will give an effective 0^ term ~ T^0^, corresponding to a stabilizing 
interaction with n = 2. In both of these interactions we have effectively massless fields 
coupling to 0, so that the time dependence of the stabilizing interaction is due to the 
red-shifting of the energy of the thermal fiuctuations of the light fields. An alternative 
is to consider the case where x is massive and is coherently oscillating in an effectively 

potential about x = 0. We assume that the x-matter density is negligible compared 
with the conventional dark matter density. In this case the x oscillation amplitude 
will be proportional to a~^^'^, resulting in a stabilizing interaction with n = 3. These 
are just a few examples; one could also consider, for example, x to be a slowly rolling 
field which triggers the phase transition, analogous to what happens in hybrid infiation 



models |E2 



There has recently been some discussion of the possibility of distinguishing between 



models with different time-dependent equations of state uj{z). In |Q it was suggested 
that in order to distinguish between models with uj{z) and u = constant we would 
require an accuracy of less than 1% in determining the luminosity distance as a function 
of red-shift, diiz), whereas 1% is regarded as an optimistic estimate of the accuracy 
of future experiements However, in |]2^, it is suggested that with an appropriate 
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fit to Lj{z), Lj{z) = uJo + ijj\z + it is possible to distinguish between different models 
using the datasets expected from the SNAP satelhte |25|. In particular, they show 
that a toy model with cjq = —0.6 and lo\ = —0.8 can be clearly distinguished from 
constant uj models. Comparing with Figure 1, we find that for n = 2, 3, 4 the expansion 
parameters {uo,uJi) are given by (-0.6,-0.33), (-0.6,-0.70) and (-0.6,-1.05) respectively. 
So at least for n = 3 and n = 4 the VRPT scenario should be clearly distinguishible 
by SNAP, and possibly for n = 2 also, although this is not directly apparent from the 
results of |^ and requires further analysis. 

In conclusion, we have introduced an alternative model for a dynamical cosmolog- 
ical constant, based on the idea that a scalar field underwent a phase transition from 
a metastable phase at a very recent epoch, z < 1.2. This VRPT scenario results in a 
characteristic evolution of the equation of state and energy density which is quite dis- 
tinct from the case of models based on slowly rolling fields, with the pressure rapidly 
rising from = —1 and rapidly decreasing at very recent times. The solutions for 
u!^{z) and P(f,{z) are uniquely determined by the present value of cj^ and the form of the 
time dependent mass squared term in the potential. Following the phase transition, 
the energy density will consist of a time dependent vacuum energy and a negative 
pressure condensate which typically fragments to a smooth pressureless density of 
non-topological solitons, '0-axitons'. 

As with most other dynamical A models, it is implicitly assumed that the reason 
for the recent dominance of the energy density and the recent occurance of the phase 
transition is connected with anthropic selection. In general, dynamical A models re- 
quire two conditions; that the energy density in the scalar field has recently become 
dominant and that the scalar field energy density is varying significantly on time scales 
of the order of H^^. Therefore two tunings are generally required. In this regard the 
VRPT scenario is on the same footing as other dynamical A models. Ultimately the 
question of which dynamical A model is correct is a matter to be decided by obser- 
vations. A future goal will therefore be to understand the detailed predictions of this 
class of dynamical A model for observable quantities; high-z supernova, quasar lensing 
statistics, the galaxy clustering power spectrum and in particular the angular CMB 
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spectrum. We expect that these will be clearly distinguishable from other dynamical 
A models, and indeed we have shown by comparison with recent analyses that at least 
some VRPT models arc likely to be distinguishable by SNAP. In addition, the idea 
of a very recent phase transition leads to other interesting issues. One is whether the 
0-axiton density could be observationally or experimentally distinguished from con- 
ventional cold dark matter. This will depend on the (p mass and on how strongly the 
(f) field interacts with ordinary matter. In addition, it is possible that the VRPT could 
result in the very recent formation of domain walls or global strings, depending on the 
symmetry associated with the potential and whether (p is real or complex. It would 
be interesting to consider whether there are any observable effects associated with the 
very recent creation of topological defects. With the possibihty of recently formed 
topological and non-top ological solitons, the phenomenology of the VRPT scenario 
for the dynamical cosmological constant may be surprisingly rich. 
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